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Abstract 

Nearest-neighbor-interaction Ising spin glasses are studied on three different hierar- 
chical lattices, all of them belonging to the Wheatstone-Bridge family. It is shown 
that the spin-glass lower critical dimension in these lattices should be greater than 
2.32 . Finite-temperature spin-glass phases are found for a lattice of fractal dimen- 
sion D ~ 3.58 (whose unit cell is obtained from a simple construction of a part of 
the cubic lattice), as well as for a lattice of fractal dimension close to five. In the for- 
mer case, the estimates of spin-glass critical temperatures associated with symmetric 
Gaussian and bimodal distributions are very close to recent results from extensive nu- 
merical simulations carried on a cubic lattice, suggesting that whole phase diagrams 
presented, obtained for couplings following non-centered distributions ~ not known 
up to the moment for Bravais lattices - should represent good approximations. 
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The Ising spin glass (SG) [l, 2| represents a very controversial model in the litera- 
ture of magnetic systems and has called the attention of many workers along the last 
decades. Theoretically, its simple formulation in terms of binary variables has led to a 
large amount of studies from both computational and analytical points of view. The 
analytical approaches have been mostly applied to its mean-field formulation, defined in 
terms of infinite-range interactions. According to its mean-field solution, the SG phase 
y described by an infinite number of order parameters (i.e., an order-parameter 



IS proper 
function 



3|]), manifesting the property of replica-symmetry breaking ^|. A major question 
in the Ising SG problem nowadays concerns the identification of which properties from 
this mean- field solution persist in the limit of short-range (e.g., nearest-neighbor) inter- 
actions. In this later case, the majority of works were concentrated on three-dimensional 
Ising SG models, for which it is generally accepted nowadays that a SG phase occurs at 



finite temperatures |4j-ll3| ; it should be mentioned that some of these studies have required 
extensive computational efforts. 

Hierarchical lattices were introduced within the context of the real-space 
renormalization-group (RG) approach, bringing the great advantage that such technique 



becomes exact for pure systems defined on these lattices These lattices are con- 

structed by carrying successive similar operations at each hierarchical level, e.g., at each 
level one replaces bonds by well-defined unit cells; typical examples of unit cells are 
presented in Fig. [H all of them belonging to the Wheatstone-Bridge (WB) family of hi- 
erarchical lattices. Two important characteristics of a given hierarchical lattice are its 
fractal dimension D and its scaling factor 6, which is defined as the smallest distance 
(counted in bonds) between the two external sites of its basic cell. The WB hierarchical 
lattices are generated by starting the process from the 0th level of the lattice-generation 
hierarchy with a single bond joining the external sites (denoted by /i and z/); then, in 
each iteration step one replaces a bond by its corresponding unit cell, in such a way that 
in its first hierarchy the lattice is represented by this unit cell. This procedure is shown 
explicitly in Fig. O for the first three levels of the WB hierarchical lattice with fractal 
dimension D = (In 5)/ (In 2) ^ 2.32. The process is continued with the hierarchical lattice 
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being constructed up to a given A/'-th hierarchy (A/" ^ 1). 

Although the RG may not be considered in general as an exact procedure for random 
systems on hierarchical lattices, it is expected to represent a good approximation, since 
in many cases, pure systems appear as particular limits of some random models. In,tlie 
particular case of Ising SGs, the hierarchical lattices have been a very useful tool 



13 



15|-|28|. essentially due to the possibility of performing relatively low-time-consuming 



numerical computations. A significant gain in the computational time in the study of SGs 
represents the main advantage of the hierarchical lattices with respect to Monte Carlo 
simulations performed on Bravais lattices. In fact, one may obtain whole phase diagrams 
for non-symmetric distributions of the couplings (e.g., by changing the probability p of 
a bimodal distribution) within typically, a few computation hours of a standard desk 
computer; on the other hand, the criticality of a given particular case (e.g., the study 
of the case p = 1/2 of the bimodal distribution for the couplings) through Monte Carlo 
simulations may require several days on high-performance computer networks. 

The Migdal-Kadanoff (MK) family of hierarchical lattices 29|, |30| have been widely 
explored, essentially due to the fact that its basic cell is formed by one-dimensional parallel 
paths each with a given scaling factor 6; its fractal dimension may be varied either by 
changing the scaling factor, or by a simple operation of adding more paths to the cell. As a 



consequence, most o 
lattices js, 8, 
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the hierarchical-lattice SG studies so far were concentrated on these 
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24| . In spite of its simplicity, some of the results obtained 



in these lattices are quite impressive: (i) The bounds for the lower critical dimension, 
which is accepted nowadays to be greater than 2, but smaller than 3, were obtained 



si about a decade before their confirmation through distinct numerical 
-|9|. In fact, the lower critical dimension on MK lattices was estimated 



on MK lattices 
investigations [6 

to be very close to 2.50 [17[; (ii) The SG critical-temperatures on the MK hierarchical 
attice of fractal dimension D = 3, for symmetric Gaussian and bimodal distributions 
st, present relative discrepancies of about 4%, taking into account the error bars, when 

compared with the recent estimates from Monte Carlo simulations on a cubic lattice 11 |. 
However, the MK lattices represent the simplest types of hierarchical lattices and 

some of its results may be quantitatively poor approximations to well-known estimates 

on Bravais lattices. Other hierarchical lattices, which present connections between such 
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parallel paths, may represent better approximations and in some cases, one may obtain 
precise estimates when compared to those known for Bravais lattices; this motivates the 
study of SG models on such lattices. In spite of this, only a few works have considered 
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25-m 



Ising SGs on hierarchical lattices different from those of the MK family 
and some of them have yielded important and stimulating results: (i) Studies on a special 
hierarchical lattice with fractal dimension D = 2 led to an estimate for the stiffness 
exponent y 19(| {y = —l/u, where u is the exponent associated with the divergence of the 



correlation length at zero temperature) in agreement with those obtained from other, more 
time-consuming, numerical approaches on a square lattice. An analysis of the ±J Ising 
SG model 25| on the same hierarchical lattice gave a ferromagnetic-paramagnetic critical 



frontier that should be a good approximation for the one of the corresponding model on a 
square lattice; (ii) Estimates of locations of multicritical points for SGs on some non-MK 
hierarchical lattices are in very good agreement with other numerical results on Bravais 



25 



27| 



lattices, as well as with a conjecture based on gauge theory 

In the present work we investigate Ising SGs defined on WB hierarchical lattices 
through the Hamiltonian, 



n 



<ij> 



The lattices generated by the cells shown in Figs, [^a), (b), and (c) present fractal di- 
mensions D = 2.32, 3.58, and 4.81; one should remind that these cells may be obtained 
from simple constructions of parts of Bravais lattices, namely, the square, cubic and four- 
dimensional hypercubic lattices, respectively jl4|. The {Jij} denote random couplings 
between two spins located at nearest-neighboring sites i and j of this hierarchical lattice 
and may follow either the Gaussian, or bimodal (±J) distributions. 



V27fj2 



exp 



2J2 



(2) 



P( J,,) = p SiJ.j - J) + (1 - 5( J,, + J) . 



(3) 
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Herein, we will restrict ourselves to those parts of the phase diagrams associated with 
ferromagnetic and SG orderings only, i.e., Jq > in Eq. ([2]),orl/2<p<lin Eq. (|3]). 

The RG procedure works in the inverse way of the lattice generation, i.e., through a dec- 
imation of the internal sites of a given cell, leading to renormalized quantities associated 
with the external sites. Defining the dimensionless couplings, Kij = f3Jij [f3 = l/{kBT)], 
the corresponding RG equations may be written in the general form. 



where Zs^^s^ represent partition functions associated with the Hamiltonian 1-i for a given 
unit cell with the external spins kept fixed {S^, = ±1), 



The RG scheme is carried by following numerically the probability distribution associated 
with the the dimensionless couphngs {Ky} [5]. Operationally, this probability distribution 
is represented by a pool of M real numbers (M is kept fixed throughout the whole RG 
procedure), from which one may compute its associated moments, at each renormalization 
step; in the limit M — )■ oo these moments should approach those of the distribution 
associated with {Kij}. The process starts by creating a pool with M coupling constants 
{Jij} generated according to one of the distributions of Eqs. (E]) or ([3]), yielding an initial 
pool of dimensionless couplings, {Kij} = f3{Jij}, for a given temperature. An iteration 
consists in M operations, where in each of them one picks randomly a set of numbers 
from the pool (each number is assigned to a bond for a given cell in Fig. [T]) in order to 
generate the effective coupling according to Eq. which will correspond to an element 
of the new pool. Following this procedure, one gets a new pool with the same size M of 
the previous one, representing the renormalized probability distribution. During the RG 
procedure, the average, < Kij >, and the width, ax =< {Kij— < Kij >)^ >^/^, are of 
particular interest for the identification of the phases, in such a way that one may obtain 
the Paramagnetic (P), Ferromagnetic (F), and Spin-Glass (SG) phases, as dominated by 
the attractors. 




(4) 



Zs^,s. = Tr{5^ (i^^,^)} [exp(-/3'H)] . 



(5) 
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< Kij > — 7- ; (Ta' — 7- ; P phase , 

< Kij > — > oo ; (Tk ^ CO (< Kij > /ctk ^ oo) ; F phase , (6) 

< Kij > — 7- ; ax ^ oo ] SG phase . 

Strictly speaking, this procedure should be carried for many different initial pools of real 
numbers, over which one may compute sample averages. However, one may also get good 
critical-frontier estimates by analyzing a sufficiently large single pool; the results that 
follow were obtained by considering a single pool of size M = 400000 real numbers. 

As expected, we did not find a finite-temperature SG phase for the hierarchical lattice 
defined by the unit cell of Fig. [D^a), for either one of the distributions of Eqs. ([2]) or 
(|3]). The phase diagram for the case of the Gaussian distribution for the couplings is 
exhibited in Fig. [3l where only two phases are present, namely, the P and F phases. 
The points used in order to draw this phase diagram were calculated from the standard 
narrowing RG procedure, with at least a two-decimal-digit certainty (error bars on third 
decimals). Due to the duality property of this unit cell one should obtain the exact 
critical temperature of the ferromagnetic Ising model on the square lattice from the phase 
diagram of Fig. [3l indeed, computing the slope of the critical frontier of Fig. [3] for {Jo/ J) 
large, one gets {k-QTc/ Jq) = 2.270 ± 0.002. Although this model does not exhibit a SG 
phase at finite temperatures, one may still compute the stiffness exponent y, which rules 
the zero-temperature scaling behavior of the width of a continuous coupling distribution 

n 

associated with blocks of linear size L [8[ , 

f{L) ^JL^- L = 2^ ■ J=< {Jij- < J,j >)2 >i/2 . (7) 

For sufficiently small values of [Jq/J), the sign of the stiffness exponent y is directly 
associated with the low-temperature phase; for a positive (negative) y the system scales to 
strong (weak) couplings, characteristic of a SG (paramagnetic) state at low temperatures. 
Therefore, one expects y < for the WB hierachical lattice of Fig. [TJ^a), whereas one 
should get y > for those of Figs. [U^b) and^c). In the first case one may use the scaling 
relation, u = —1/y [8], for a phase transition in the limit T — )■ 0, where u is the exponent 
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associated with the divergence of the correlation length, ,^ ~ T ^. We have computed 
y = -0.290 ± 0.003, leading to u = 3.45 ± 0.04, for all Jq in the interval < (Jq/J) < 1 
(cf. Fig. [3]). This estimate coincides, within the error bars, with those obtained for the 



particular case Jq = in Ref. 



19| for the same hierarchical lattice, as well as with several 



estimates for the square 
and y = -0.284 ± 0.004 



attice, like y = -0.291 ± 0.002 [8], y = -0.287 ± 0.004 |3l|. 



32|. Although the lattice considered herein presents a fractal 



dimension D ^ 2.32, which is about 15% higher than the dimension of the square lattice, 
it has been much used in the literature as an approach for models on the later, essentially 
because it may be obtained through a simple construction from a piec e of the square 
lattice, keeping one of its most important properties, the self-duality [14 1. 

The phase diagram for the case of the hierarchical lattice defined by the unit cell of 
Fig. lU^a), with a bimodal distribution for the couplings, is qualitatively similar to the 
one already investigated in Ref. 25(|. From the present analysis one concludes that the 
SG lower critical dimension in the WB hierarchical-lattice family should be greater than 
2.32; this result is in agreement with calculations carried on MK lattices, where this lower 
critical dimension was estimated to be very close to 2.50 [l^. 

The phase diagrams for the hierarchical lattice defined by the unit cell of Fig. Wljo) 
are shown in Fig. Hj for couplings following the Gaussian distribution [Fig. HJ^a)], or 
the bimodal distribution [Fig. Hl^b)]. The coordinates of the most important critical 
points of these phase diagrams are given explicitly in Tables I and II. Similarly to the 
previous case, the points used in order to draw these phase diagrams were calculated 
from the standard narrowing RG procedure, with at least a two-decimal-digit certainty 
(error bars on third decimals). It is important to notice the estimates of the critical 
temperatures for symmetric distributions, namely, (ksTc/J) = 0.980(2) and (kBTc/J) = 
1.112(2); these estimates are in very good agreement with recent results from Monte 
Carlo simulations on a cubic lattice, which yielded {kBTd J) = 0.951(9) and (ksTc/J) = 
1.120(4), for symmetric Gaussian and bimodal distributions, respectively [11]. Comparing 
the present results with those of Ref. Uj , taking into account the error bars, one finds a 
relative discrepancy of about 2% in the Gaussian case, whereas for the symmetric bimodal 
distribution the two estimates essentially coincide (leading to a relative discrepancy of 
about 0.2%). Such good agreements on critical-temperature estimates, between those of 
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a hierarchical lattice with fractal dimension D ^ 3.58 and those of the cubic lattice, are 
justified by the fact that the cell of Fig. [I](b) may be obtained from a simple construction 
of a piece of the cubic lattice [l^. On the other hand, the critical temperature associated 
with the corresponding ferromagnetic Ising model, which may be obtained either from the 
p = 1 critical point separating phases P and F in Fig. IDj^b), or by computing the slope 
of the critical frontier of Fig. Hl^b) for {Jo/ J) large, is given by {k-oTc/ Jq) = 5.46 ± 0.01, 
which is somewhere in between the estimates of critical temperatures for the ferromagnetic 
Ising model on the cubic [{k-QTc/ Jq) ~ 4.512] and the four-dimensional hypercubic lattice 
[(/cbTc/Jq) ~ 6.680], obtained from Monte Carlo simulations [33|. 

Another important aspect of the phase diagrams shown in Fig. H] is a small reentrance 
effect in the region slightly to the right of the multicritical point; by lowering the temper- 
ature, one may go from a high-temperature disordered (P) phase to an ordered one (F 
phase) and then, back to a more disordered (SG) phase at low temperatures. By compar- 
ing the coordinates of the multicritical points and the respective zero-temperature ones in 
Tables I and II, one sees that this effect, although very weak, is outside the error bars of 
the method. This is illustrated in Fig. |5l which shows the evolution of the average < Kij > 
with the RG step n, in the case of couplings following an initial Gaussian distribution, 
for typical values of Jq/J around the critical frontier separating the phases SG and F. 
In Fig. [5t^a) one sees the behavior of < Kij > for a fixed temperature [{kBT/ J) = 0.95], 
below the SG critical temperature, showing that for {Jq/J) = 0.556 one is still on the 
SG phase, whereas < K^j > diverges for {Jq/J) = 0.557, signalling a F phase. The 
same procedure is applied to zero temperature in Fig. |5]^b), where one sees these two 
different types of behavior for {Jq/J) = 0.566 and {Jq/J) = 0.567, respectively. It should 
be mentioned that this type of reentrance phenomena has been observed experimentally. 



e.g., in EuxSri_xS |3J] and AuFe (14% Fe) [35|. Theoretically, it has been also found in 
the mean-field replica-symmetric solution (which turns out to be unstable at low temper- 
atures) of the Ising SG, although its correct solution, characterized by replica-symmetry 
breaking, has washed away the reentrance jl-3]. It is important to stress that reentrance 
phenomena on theoretical SG models are very subtle and difficult to be obtained numer- 



ically, in such a way that on. 
effect so far (see, e.g., Refs. 



y a few works in the literature were able to capture such an 
25|,|36|). 
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Hierarchical 
lattice 


kBTc/J 
(Jo = 0) 


Jo/J 
(T = 0) 


Multicritical Point 


Cell 1(b) 
Cell 1(c) 


0.980(2) 
2.35(1) 


0.5665(5) 
0.374(3) 


{kBTjJ) = 1.690(2); (Jq/J) = 0.538(2) 
{kBTjJ) = 2.77(1); (Jo/J) = 0.35(1) 



TABLE I: Values of typical critical points in the phase diagrams of the Ising SG, on the hier- 
archical lattices defined by unit cells of Figs. [Hb) andlljc), for couplings following a Gaussian 
distribution [Eq. ([2])]: critical temperature for a symmetric distribution, zero-temperature criti- 
cal point separating phases SG and F, and coordinates of the multicritical point where the three 
critical frontiers meet. The error bars refer to the usual approach to criticality characteristic of 
the RG technique. 

We have also computed the zero-temperature stiffness exponent y for the phase diagram 
of Fig. m^a). As expected, a positive value, y = 0.297 ± 0.003, signals a SG phase at 
finite temperatures; furthermore, this estimate is universal for all values of Jq within the 
SG phase. This value is close to the recent one obtained for a MK lattice of fractal 
dimension Z) = 3 (y ~ 0.27 [13]), but it is slightly larger than those for the cubic lattice, 
y = 0.19 ± 0.01 [8] and y = 0.20 ± 0.05 Q- 

The phase diagrams for the hierarchical lattice defined by the unit cell of Fig. [U^c), 
for couplings following either initial Gaussian of bimodal distributions, are qualitatively 
similar to those shown in Figs. Hl^a) andlD^b), respectively. Values of typical critical points 
in these phase diagrams are given explicitly in Tables I and II. The higher complexity of 
this lattice generates larger errors in the numerical estimates; however, the points in such 
phase diagrams were computed with at least a one-decimal-digit certainty (error bars on 
second decimals). In particular, one should mention the estimates of the critical temper- 
atures for symmetric distributions, namely, (fc^Tc/J) = 2.35(1) and (fc^Tc/J) = 2.515(2), 
which are about 25% above those obtained from Monte Carlo simulations for the four- 
dimensional hypercubic lattice, i.e., {kBTc/ J) ~ 1-80, for the Gaussian distribution js^], 



and (ksTc/J) ~ 2.00, for the bimodal distribution |38|, respectively. These overestimates 
are consistent with the fact that the fractal dimension of the lattice considered herein is 
well above four {D ^ 4.81). Although we are not aware of numerical simulations carried 
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Hierarchical 
lattice 


ksTc/J 
{p=l/2) 


Pc 

(r = o) 


Multicritical Point 


Cell 1(b) 
Cell 1(c) 


1.112(2) 
2.515(2) 


0.760(1) 
0.667(2) 


{kBTjJ) = 1.620(2); = 0.745(2) 
(kBTc/J) = 2.836(2); = 0.664(2) 



TABLE II: Values of typical critical points in the phase diagrams of the Ising SG, on the hier- 
archical lattices defined by unit cells of Figs. [Hb) andfljc), for couplings following a bimodal 
distribution [Eq. critical temperature for a symmetric distribution, zero-temperature criti- 
cal point separating phases SG and F, and coordinates of the multicritical point where the three 
critical frontiers meet. The error bars refer to the usual approach to criticality characteristic of 
the RG technique. 



for Ising SGs on lattices of higher dimensions, the present estimate for the symmetric 
bimodal distribution is slightly below the one obtained through series expansions on the 
five-dimensional hypercubic lattice, (fc^Tc/J) ^ 2.57 jsQ]. 

To conclude, we have investigated nearest-neighbor-interaction Ising spin glasses de- 
fined on three different hierarchical lattices, belonging to a family of Wheatstone-Bridge 
lattices, characterized by fractal dimensions D ^ 2.32, 3.58, and 4.81. The interactions 
among pairs of spins were chosen from either a Gaussian, or a bimodal (±J), non-centered 
distribution. Through calculations of the stiffness exponent, we have shown that the 
spin-glass lower critical dimension in these lattices should be greater than 2.32 . Finite- 
temperature spin-glass phases were found for the lattices of fractal dimension D ^ 3.58 
and D ^ 4.81. In the former case, whose hierarchical lattice is constructed as an approx- 
imation to the cubic lattice, the estimates of spin-glass critical temperatures associated 
with zero-centered Gaussian and bimodal distributions are very close to the most recent 
estimates from extensive numerical simulations carried on the cubic lattice. Phase dia- 
grams were obtained for couplings following non-centered distributions - which are not 
entirely known, up to the moment, for spin-glass models on Bravais lattices - either for 
Jo > in the case of the Gaussian distribution, orl/2<j9<lin the case of the bimodal 
distribution; these phase diagrams are expected to represent good approximations for the 
corresponding models on Bravais lattices. 
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Figure Captions 

Fig. 1: Basic cells of the WB family of hierarchical lattices, with a scaling factor b = 2. 
In each case, the empty circles (/x and u) represent the external sites of the cell, whereas 
the black circles are internal sites to be decimated in the RG procedure, (a) The WB 
cell of fractal dimension D = (In5)/(ln2) ^ 2.32; (b) The WB cell of fractal dimension 
D = (In 12) /(In 2) ^ 3.58; (c) The WB cell of fractal dimension D = (In 28) /(In 2) ^ 4.81. 
For clearness, all the bonds of cell (c) are not shown explicitly; the external site /i is 
connected directly to sites 5, 6, 7, and 8 (as illustrated by the bond connecting sites /i and 
5, whereas site u is connected directly to sites 1, 2, 3, and 4 (as illustrated by the bond 
connecting sites u and 3). 

Fig. 2: First three levels (levels 0, 1, and 2) in the generation process of the Wheatstone- 
Bridge hierarchical lattice with fractal dimension D ^ 2.32. The lattice-generation process 
starts at the 0-th level with a single bond and at each step a bond is replaced by its 
respective basic cell. In each level, the empty circles {fi and u) represent the external sites 
of the lattice, whereas the black circles are internal sites to be decimated in successive 
RG operations. 

Fig. 3: Phase diagram for an Ising SG on the hierarchical lattice defined by the unit cell 
of Fig. [U^a), with a Gaussian distribution for the couplings. 

Fig. 4: Phase diagrams of the Ising SG on the hierarchical lattice defined by the unit 
cell of Fig. IH^b): (a) Gaussian distribution for the couphngs; (b) bimodal distribution for 
the couplings. 

Fig. 5: Evolution of the average of couplings with the RG step n, for the case of 
an initial Gaussian distribution, considering typical values of Jq/J around the critical 
frontier separating the phases SG and F. (a) Fixed temperature below the SG critical 
temperature; (b) Zero temperature. 
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FIG. 1: Basic cells of the WB family of hierarchical lattices, with a scaling factor 6 = 2. In each 
case, the empty circles (/i and u) represent the external sites of the cell, whereas the black circles 
are internal sites to be decimated in the RG procedure, (a) The WB cell of fractal dimension 
D = (In5)/(ln2) « 2.32; (b) The WB cell of fractal dimension D = (Inl2)/(ln2) ^ 3.58; (c) 
The WB cell of fractal dimension D = (In 28)/ (In 2) ^ 4.81. For clearness, all the bonds of cell 
(c) are not shown explicitly; the external site /x is connected directly to sites 5, 6, 7, and 8 (as 
illustrated by the bond connecting sites fi and 5, whereas site u is connected directly to sites 
1,2,3, and 4 (as illustrated by the bond connecting sites u and 3). 
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FIG. 2: First three levels (levels 0, 1, and 2) in the generation process of the Wheatstone-Bridge 
hierarchical lattice with fractal dimension D ~ 2.32. The lattice-generation process starts at 
the 0-th level with a single bond and at each step a bond is replaced by its respective basic cell. 
In each level, the empty circles (/x and u) represent the external sites of the lattice, whereas the 
black circles are internal sites to be decimated in successive RG operations. 
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FIG. 3: Phase diagram for an Ising SG on the hierarchical lattice defined by the unit cell of 
Fig. dl^a), with a Gaussian distribution for the couplings. 




FIG. 4: Phase diagrams of the Ising SG on the hierarchical lattice defined by the unit cell of 
Fig. [T]||b): (a) Gaussian distribution for the couplings; (b) bimodal distribution for the couplings. 
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FIG. 5: Evolution of the average of couplings with the RG step n, for the case of an initial 

Gaussian distribution, considering typical values of Jq/J around the critical frontier separating 
the phases SG and F. (a) Fixed temperature below the SG critical temperature; (b) Zero 
temperature. 



